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The  concept  of  adaptive  filters  and  adaptive 


communication  systems  is  destined  to  play  a  major  role  in 
the  future.  With  this  thought  in  mind,  I  embarked  on  a 
study  that  analyzed  and  simulated  an  adaptive  matched  filter 
that  uses  the  Griffiths  algorithm.  The  original  idea  came 
from  readings  offered  to  me  by  Martin  DeSimio  of  Wright- 
Patterson  Air  Force  Base,  Ohio.  This  idea  blossomed  into 
this  thesis,  which  is  a  result  of  approximately  nine  months 
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This  thesis  presents  a  CW  and  noise  jamming 
analysis  of  an  adaptive  matched  filter  that  (1)  uses  the 
Griffiths  algorithm  and  (2)  has  a  pseudonoise  sequence  as  an 
input.  The  analysis  is  conducted  over  several  jamming 
powers,  frequencies,  and  phases.  The  Griffiths  adaptive 
matched  filter  is  shown  to  converge  for  raised  cosine  pulses 
that  experience  no  distortion,  quadratic  delay  distortion, 
and  cubic  delay  distortion.  The  Griffiths  adaptive  matched 
filter  diverges  for  pulses  that  experience  linear  delay 
distortion  even  though  the  convergence  rate  constant  is 
within  limits.  Throughout  the  analysis  the  P-vector  is 
determined  apriori  and  held  constant.  The  Griffiths  filter 
is  shown  to  converge  for  CW  jamming  and  noise  jamming. 

Noise  jamming  is  shown  to  be  more  effective  in  the  higher 
power  ranges.  A  comparison  is  made  between  the  Griffiths 
adaptive  matched  filter  and  an  adaptive  matched  filter  that 
uses  the  LMS  algorithm.  The  degradation  in  performance  of 
the  Griffiths  filter  compared  to  an  LMS  filter  that  uses  a 
stored  reference  is  calculated  for  several  selected  runs. 

The  actual  computer  programs  used  are  presented  in  the 
Appendix. 
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PERFORMAKCE  OF  AN  ADAPTIVE 
MATCHED  FILTER  USING  THE 
GRIFFITHS  ALGORITHM 
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CHAPTER  1  :  PRELIMINARIES 


1-1  INTRODUCTION 

This  chapter  addresses  basic  preliminaries.  A 
justification  for  this  research  effort  is  stated  along  with 
selected  background  information.  The  concept  of  adaptive 
matched  filtering  is  explained  along  with  differences  from 
several  related  concepts.  A  statement  of  the  problem  is 
presented,  and  the  scope  of  the  research  effort  is  defined. 


1-2  JUSTIFICATION 

Performance  of  a  military,  digital  communications 
receiver  in  a  jamming  or  intentional  interference 
environment  is  a  primary  consideration  for  military 
planners.  Unexpected  performance  of  a  receiver  on  the 
electronic  battlefield  can  have  disastrous  consequences  for 
the  user.  The  Air  Force  has  realized  the  importance  of 


I 

I 

I 

. 
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secure  and  effective  electronic  communications  and  has 


incorporated  the  concept  into  its  basic  aerospace  doctrine. 


Commanders  rely  on. .. secure. . . 
communications.  Communications  are 
the  means  through  which  a  commander 
transmits  and  receives  information 
about  the  enemy,  coordinates  with 
friendly  forces,  and  commands  and 
controls  assigned  forces  (1:2-21). 


The  criticality  of  good 
communications  to  military  operations 
cannot  be  overemphasized.  Secure,  jam 
resistant  communications  must  be 
established  (2:8). 


Since  secure  and  jam-resistant  communications  are 
important  to  the  Air  Force,  various  electronic  communication 
receivers  have  been  deployed  to  counteract  the  expected 
jamming  threat.  These  receivers  use  many  different  methods; 
however,  one  trend  has  been  toward  adaptive  signal 
processing.  This  trend  toward  adaptive  signal  processing  is 
expected  since  the  electronic  battlefield  is  a  changing 
environment,  and  one  can  not  predict  what  type  of  electronic 
signals  will  be  present. 


Adaptivity  is  needed. . .because  the 
environment  is  apriori  unknown,  time 
varying,  and  noncooperative.  The  present 
day  systems  have  fixed  architecture  and  a 
limited  degree  of  adaptivity.  This 
corresponds  to  a  relatively  simple  hardware 
but  provides  non-optimal  performance  when 
the  environment  differs  from  that  assumed 
in  the  radar  design. .. .Adaptivity  should 
be  used  extensively  in  each  subsystem  of 
the  radar  (3:178). 
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1-3  BACKGROUND 

One  use  of  adaptive  signal  processing  to  help 
counter  interference  or  jamming  is  in  the  matched  filter 
circuitry  of  a  digital  communications  receiver.  In  this 
situation,  fixed  coefficient,  tapped-delay-line  matched 
filters  are  designed  so  that  they  can  adapt  to  received 
waveforms.  This  would  be  beneficial  since  one  could  not 
predict  what  type  of  waveform  would  be  received  in  a  jamming 
or  interference  environment.  Once  the  matched  filters  have 
adapted  to  the  received  waveform,  the  standard  decision 
circuitry  and  test  statistics  that  were  used  in  the  fixed 
coefficient  case  can  then  be  used. 

The  ideal  adaptive  matched  filter  adapts  itself  to 
the  received  waveform  so  that  the  filter  becomes  matched  to 
the  received  waveform.  The  matched  filter  is  made  to  adapt 
so  that  "the  filter  is  optimized  by  minimizing  the  mean- 
square  value  of  the  error  signal"  (10:27).  A  filter  which 
minimizes  the  mean-square  error  is  optimum  in  the  mean- 
square  sense.  The  adaptive  matched  filter  minimizes  the 
mean-square  error  between  the  actual  filter  output  and  a 
desired  signal.  By  minimizing  the  mean-square  error,  the 
adaptive  matched  filter's  output,  or  test  statistic,  can 
remain  approximately  constant  for  a  given  transmitted 
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symbol,  and  the  signal  to  noise  ratio  at  the  output  will  be 
maximized  at  the  sampling  time. 

Several  other  concepts  closely  parallel  adaptive 
matched  filtering  but  will  not  be  covered  in  this  thesis. 
These  other  concepts  are  (1)  equalization,  (2)  adaptive 
interference  suppression,  and  (3)  adaptive  beam  forming. 

In  equalization,  the  distortion  caused  by  a 
channel  is  compensated  for  so  that  intersymbol  interference 
(ISI)  can  be  reduced.  Equalization  is  used  to  combat  ISI 
(9:545).  In  equalization,  the  sampling  rate  is  the  symbol 
rate  (11:105).  By  contrast,  in  adaptive  matched  filtering 
the  sampling  rate  is  an  integer  multiple  of  the  symbol  rate. 
Adaptive  matched  filtering  is  used  to  optimally  detect  in 
the  mean-square  sense  the  received  digital  waveform,  and, 
therefore,  it  will  combat  intentional  interference  and 
noise.  It  is  not  intended  to  correct  or  eliminate  ISI. 

In  adaptive  interference  suppression,  the 
interference  is  notched  out  in  the  frequency  domain  to 
obtain  an  approximation  of  the  signal.  The  method  is  used 
in  such  devices  as  adaptive  notch  filters  and  adaptive  line 
enhancers.  The  basic  idea  is  to  suppress  the  interference 
or  noise  and  extract  the  signal  (19:1698).  By  contrast,  in 
adaptive  matched  filtering,  no  attempt  is  made  to  suppress 
the  interference  in  the  frequency  domain.  Instead,  the 
adaptive  matched  filter  maximizes  the  signal  to  noise  ratio 
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at  the  output  by  minimizing  the  mean-square  value  of  the 
error  signal. 

In  adaptive  beam  forming  an  antenna  radiation 
pattern  is  purposely  nulled  out  in  specific  directions  in 
order  to  suppress  interference  (8:409).  The  method  works 
well  if  the  signal  and  interference  can  be  separated 
spatially.  In  adaptive  matched  filtering  the  concept  of 
look-directions  has  no  physical  meaning  and  matched 
filtering  is  not  accomplished  in  the  antenna  section.  In 
adaptive  matched  filtering  the  suppression  of  interference 
is  accomplished  by  maximizing  the  signal  to  noise  ratio  at 
the  receiver  decision  circuitry. 

The  concepts  of  equalization,  adaptive 
interference  suppression,  adaptive  beam  foirming,  and 
adaptive  matched  filtering  all  have  basically  the  same  goal. 
They  all  try,  in  some  way,  to  provide  reliable  communication 
in  a  hostile  (or  less  than  optimum)  environment.  The  actual 
concept  employed  will  depend  on  the  expected  threat  and  the 
decisions  made  by  senior  commanders  (20:28).  In  this 
thesis,  the  concept  of  adaptive  matched  filtering  will  be 
covered.  Also,  the  performance  of  an  adaptive  matched 
filter  against  an  expected  jamming  threat  will  be  covered. 

An  adaptive  matched  filter  is  a  finite  impulse 
response  (FIR)  filter  whose  impulse  response  is  adjusted  in 
some  optimum  way.  The  procedure  for  adjusting  the  adaptive 
matched  filter  is  known  as  the  adaptive  algorithm.  Two  of 
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the  more  popular  algorithms  are  the  conventional  least-mean- 
square  (IKS)  algorithm  developed  by  Widrow  and  Hoff  and,  a 
IKS  variant,  the  Griffiths  algorithm.  This  thesis  is  mainly 
concerned  with  the  Griffiths  algorithm. 

The  Griffiths  algorithm,  which  is  also  known  as 
the  P-vector  algorithm  or  modified  IKS,  was  developed  by 
L.J.  Griffiths  in  1969  and  has  an  advantage  in  that  the 
desired  signal  term  is  not  used  directly  in  the  algorithm 
(4:1696).  Since  the  desired  signal  is  not  used  directly  in 
the  algorithm,  the  requirement  of  the  receiver  knowing  the 
desired  signal  (i.e  having  a  stored  reference)  is  no  longer 
necessary.  Also,  Feuer  and  Weinstein  have  shown  that  the 
Griffiths  algorithm  produces  lower  misadjustment  (see 
Appendix  B)  than  the  conventional  IKS  algorithm  when  the 
correlation  is  low  between  the  received  data  and  the  desired 
signal  (5:226) . 

Due  to  its  lower  misadjustment  compared  to  the 
conventional  IKS,  the  Griffiths  algorithm  is  showing 
somewhat  increased  usage  in  the  reception  of  signals  in  the 
low  SNR  environment.  A  recent  trend  has  been  to  use  the 
Griffiths  adaptive  algorithm  for  adaptive  processing  of 
signals  in  the  matched  filter  circuitry  of  digital 
communication  receivers. 


The  Griffiths  algorithm  is  being  implemented  in 


the  matched  filter  circuitry  of  digital  communication 
receivers  due  to  its  expected  better  misadjustment  in  a  low 
SNR  environment;  however,  no  comprehensive  study  has  been 
accomplished  concerning  the  performance  of  the  Griffiths 
algorithm  in  a  jamming  or  intentional  interference 
environment  with  direct  application  to  a  communications 
system.  The  question  that  needs  to  be  answered  is  whether 
the  Griffiths  adaptive  matched  filter  provides  any 
advantages  over  conventional  methods  (i.e.  LMS)  in  a 
communications  environment.  For  example,  Feuer  and 
Weinstein's  work  investigated  performance  in  Gaussian  noise 
only. 

Since  performance  in  a  jamming  environment  is 
important,  the  analysis  of  a  Griffiths  adaptive  filter  in  a 
jamming  environment  (with  direct  applications  to 
communications)  needs  to  be  accomplished.  This  thesis  will 
investigate  performance  against  a  continuous  wave  (CW) 
jammer  and  a  noise  jammer. 


( 
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This  thesis  will  describe  a  simulation,  via 


computer  program,  of  an  adaptive  matched  filter  that  uses 


the  complex  Griffiths  algorithm.  The  following  objectives 
will  be  accomplished  in  this  thesis: 


A.  Perform  a  theoretical  analysis  of  the 
Griffiths  adaptive  matched  filter. 

B.  Model  the  Griffiths  adaptive  matched 
filter  via  a  computer  program. 

C.  Input  a  pseudonoise  (PN)  test  signal  of 
length  31  into  the  Griffiths  adaptive  matched  filter  and 
analyze  performance.  Performance  will  be  analyzed  over  a 
range  of  jamming  to  signal  power  ratios  (-10  DB,  0  DB,  and 
10  DB)  and  signal  phase  distortion  levels.  Signal  phase 
distortion  levels  will  be  one  of  the  following  four  types: 

1)  No  distortion 

2)  Linear  delay  distortion 

3)  Quadratic  delay  distortion 

4;  Cubic  delay  distortion 

A  comparison  of  results  to  the  LMS  adaptive  matched  filter 
will  be  accomplished.  The  jammer  will  be  modeled  as  both  a 
CW  jammer  and  noise  jammer. 
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CHAPTER  2  :  FUNDAMENTAL  CONCEPTS 


2-1  INTRODUCTION 

In  this  chapter  the  fundamental  concepts  of 
Griffiths  adaptive  matched  filtering  will  be  explained. 
First,  the  simple  definition  of  a  digital  filter  will  be 
addressed  and  then  the  more  advanced  concepts  of  matched 
filtering  and  adaptive  filtering  will  be  presented. 
Finally,  an  explanation  of  the  complex  Griffiths  algorithm 
will  be  presented  along  with  its  relationship  to  the 
conventional  LMS  algorithm. 


2-2  BASIC  DIGITAL  FILTERING 

Basic  digital  filtering  is  concerned  with  the 
extraction  of  infozrmation  from  some  type  of  data.  In  a 
digital  communication  system,  the  information  is  usually 
extracted  from  a  received  data  sequence.  The  received  data 
sequence  contains  a  signal  of  interest  corrupted  by  noise. 
The  digital  filter  operates  on  the  received  data  sequence  in 
order  to  make  an  estimate  of  the  signal  of  interest.  In 
many  cases,  such  as  matched  filtering,  the  digital  filter 
produces  a  test  statistic  which  indicates  which  symbol 
waveform  was  transmitted. 
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Basic  digital  filters  are  usually  implemented  in 
(1)  filtering  operations,  (2)  smoothing  operations,  or  (3) 
prediction  operations  (10:1).  In  this  thesis,  only  the 
concept  of  filtering  will  be  treated.  Filtering  is  defined 
to  be  "the  extraction  of  information  about  a  quantity  of 
interest  at  time  t  by  using  data  measured  up  to  time  t" 
(10:1)  . 

Digital  filters  are  usually  represented  by  the 
block  diagram  given  in  Figure  2-1.  r(k),  which  is  a  sampled 

version  of  a  cyclostationary  (see  Appendix  B)  random 
process,  is  the  input  sequence.  r(k)  is  also  assumed  to 
have  been  sampled  at  or  above  the  Nyquist  rate.  h(k)  is 
the  filter  impulse  response,  and  y(k)  is  the  filter  output 
sequence.  The  representation  in  Figure  2-1  is  the  time 
domain  representation.  An  equivalent  frequency  domain 
representation  could  be  obtained  by  taking  the  appropriate 
Fourier  transforms  (7:101).  In  this  thesis,  only  time 
domain  representations  will  be  considered. 

Digital  filters  are  classified  as  either  infinite 
impulse  response  (HR)  or  finite  impulse  response  (FIR)  . 

This  thesis  will  be  concerned  with  causal  FIR  filters. 
Causality  restricts  the  filter  to  have  an  output  based  only 
on  present  and  past  inputs.  FIR  filters  are  also  known  as 
transveral  filters  or  tapped-delay-line  filters. 


10 


r(k)  - 

input 

sequence 


h(k) 


y(k) 

output 

sequence 


h(k)  is  the  filter's  impulse  response 


Figure  2-1  Block  Diagram  of  a  Digital  Filter 
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2-3  MATCHED  FILTERING 

In  this  thesis  we  are  concerned  with  digital 
filters  implemented  in  the  matched  filter  circuitry  of  a 
digital  communication  receiver.  “A  matched  filter  is  a 
linear  filter  designed  to  provide  the  maximum  signal-to- 
noise  power  ratio  at  its  output  for  a  given  transmitted 
symbol  waveform”  (11:88).  Figure  2-2  shows  a  typical  block 
diagram  of  a  matched  filter  receiver  for  binary  phase  shift 
keyed  (BPSK)  signals. 

In  Figure  2-2  the  received  data  vector  r(k)  is 
processed  in  two  matched  filters.  The  outputs  of  the 
matched  filters  are  sampled  at  time  t=T.  The  sampling  time 
T  corresponds  to  the  symbol  time  and  occurs  when  jr(k)  is 
completely  justified  within  the  matched  filter  (sampling 
time  is  also  related  to  filter  length  and  the  eye  opening  of 
the  eye  diagram) .  The  two  sampled  outputs  (or  test 
statistics)  are  then  compared,  and  a  decision  is  then  made 
on  which  symbol  was  sent.  Throughout  the  entire  process 
synchronization  is  assumed  to  be  maintained. 

In  the  development  of  the  matched  filter  (11:88), 
a  known  signal  8(k)  and  additive  white  gaussian  noise 
fi(k)  are  assumed  to  be  input  to  the  matched  filter 
circuitry.  Hence,  the  received  signal  £(k)  is  of  the  form 

I(k)  =  fi(K)  +  a(k)  (1) 
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Figure  2-2  Block  Diagram  for  a  Matched  Filter 
Receiver  for  BPSK 
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The  matched  filter  can  be  implemented  digitally  as 
a  tapped-delay-line.  A  tapped-delay-line  matched  filter  is 
shown  in  Figure  2-3.  In  Figure  2-3,  five  samples  occur  each 
symbol  duration  T.  The  weights  Wj^  thru  Wg  are 
determined  by  the  signal  vector  8(k)  and  are  fixed  and 
known  apriori.  The  output  of  the  tapped-delay-line  is 
defined  by  the  matrix  equation 

Y(k)  =  r'^(k)  w(k)  (2) 

Note  that  a  scalar  (or  test  statistic)  is  output  for  each 
vector  input.  When  Y(k)  is  sampled  at  time  t=T,  the 
symbol  duration  time,  the  signal  to  noise  ratio  is 
maximized.  Note  that  r(k)  ,  uW  i  s(J<)  »  KW  >  and  Y(k) 
can  all  be  complex  quantities. 


2-4  ADAPTIVE  MATCHED  FILTERING 

In  the  previous  two  sections,  discussion  has 
centered  on  basic  digital  filtering  and  matched  filtering. 

It  was  shown  that  a  matched  filter  could  be  implemented  as  a 
tapped-delay-line  FIR  filter.  When  the  signal  vector  s(k) 
is  known  exactly  and  the  noise  vector  n(k)  is  an  additive 
white  gaussian  noise  process,  the  matched  filter  can  be 
implemented  as  a  fixed  weight  tapped-delay-line. 
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Figure  2-3  Tapped-delay-line  Matched  Filter 


In  most  instances,  however,  the  signal  s(k)  is 
not  Jcnown  exactly.  Distortion  in  the  communication  channel 
may  distort  £(k)  in  amplitude  and  phase.  Also, 
intentional  interference  or  jamming  may  affect  s(k)  and 
contribute  to  uW  •  In  cases  where  fi(k)  is  distorted, 
the  fixed  weight  matched  filter  is  less  than  optimum  since, 
in  the  derivation  of  the  matched  filter,  5(k)  was  assumed 
to  be  known  exactly. 

One  way  of  overcoming  the  problem  of  nonoptimum 
performance  is  to  make  the  fixed-weight  matched  filter  adapt 
to  the  received  vector  £(k) .  This  can  be  accomplished  by 
modifying  the  diagram  of  Figure  2-3  to  that  of  Figure  2-4. 

Figure  2-4  shows  the  classic  tapped-delay-line 
modified  to  perform  the  adaptive  filtering  operation 
(10:91).  In  this  case,  the  weight  vector  w(k)  is  adjusted 
to  match  the  received  vector  £(k) .  Once  the  weights  have 
adjusted  to  r(k) ,  the  filter  becomes  matched.  When  the 
weight  vector  w(k)  has  completely  adapted,  Y(k)  will 
approximately  equal  D(k)  and  the  mean-square  value  of  the 
error  signal  will  be  minimized. 

The  advantage  of  having  an  adaptive  matched  filter 
is  that,  if  X(k)  changes  due  to  amplitude,  phase,  or 
interference  distortion,  the  filter  can  adapt  to  the  new 
X(Jc) . 
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Figure  2-4  Adaptive  Matched  Filter 
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The  adaptive  algorithm  in  Figure  2-4  is  the 
procedure  by  which  the  weights  are  updated.  In  this  thesis, 
the  complex  Griffiths  algorithm  will  be  implemented  as  the 
adaptive  algorithm. 

The  complex  Griffiths  algorithm  can  be  derived 
from  the  Widrow-Hoff  IMS  algorithm  (8:412).  The  complex 
Widrow-Hoff  LMS  algorithm  is  given  as  (15:719) 

w(k+l)  =  w(k)  +  2n  e(k)  £*(k)  (3) 

where 

k  *  the  k^^  iteration  or  clock  cycle 
S?(k+1)  =  the  updated  complex  weight  vector 
w(k)  =  the  previous  complex  weight  vector 
^  =  convergence  rate  constant 
e(k)  =  complex  error  signal 
£*(k)  =  the  complex  conjugate  of  the  received  vector 

If  the  substitution  e(k)  =  D(k)  -  Y(k)  ,  where  D(k)  is  the 
desired  signal  and  Y(k)  is  the  actual  output,  is  made  in 
equation  (3)  then 

K(k+1)  =  w(k)  +  2^1  [  D(k)  -  Y(k)  ]  r*(k) 


(4) 


After  expanding,  equation  (4)  becomes 


lC(k+l)  =  !r(k)  +  2u  D(k)  i*(k)  -  2n  Y(k)  r*(k)  (5) 

Equation  (5)  is  an  expanded  form  of  the  LMS  algorithm.  In 
1969  L.J.  Griffiths  investigated  the  possibility  of  removing 
the  D(k)  term  from  the  algorithm.  The  motivation  for 
this  was  that  in  certain  adaptive  receiving  antenna  arrays  a 
desired  signal  is  not  available  for  use  in  the  adaptive 
algorithm,  and  methods  such  as  the  LMS  algorithm,  which 
require  a  desired  signal,  cannot  be  directly  applied 
(6:705).  Griffiths'  work  led  him  to  replace  the  D(k)  £*(k) 
term  with  its  expected  value  (4:1699).  With  this 
substitution,  equation  (5)  becomes 

!?(k+l)  =  w(k)  +  2u  E{  D(k)  r*(k)  )  -  2u  Y(k)  £*(k)  (6) 

and  let  P  =  E{  D(k)  i*(k)  ) 

ir(k+l)  =  w(k)  +  2^i  £  -  2»i  Y(k)  £*(k)  (7) 

Equation  (7) ,  as  shown  above,  is  the  complex 
Griffiths  algorithm  (6:705).  The  middle  term  in  equation 
(7),  also  known  as  the  P-vector,  is  a  constant  determined 
from  apriori  information  (4:1699).  The  P-vector  in  this 
thesis  will  be  determined  from  a  signal  that  has  no 
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distortion.  The  degradation  in  filter  performance  will  be 
analyzed  when  the  actual  received  signal  has  phase 
distortion  added  (which  is  different  from  what  the  P-vector 
was  determined  from) . 

Equation  (7)  can  be  shown  to  converge  to  the  best 
mean-square  estimate  of  the  optimum  Wiener  solution  as  the 
number  of  adaptations  approaches  infinity  (8:413).  This 
convergence  is  guaranteed  as  long  as  the  convergence  rate 
constant  satisfies 

0  <  |j,  <  -i-  (8) 

’^max 

where 

Tijjax  “  largest  eigenvalue  of  the  input 
correlation  matrix 

"Like  the  LMS  algorithm,  the  Griffiths  algorithm  is  unbiased 
and  produces  at  convergence  an  expected  steady-state 
solution  that  is  the  true  least-squares  solution"  (8:413). 
The  advantages  of  the  Griffiths  algorithm  are  that  the 
desired  signal,  D(k) ,  is  not  directly  used  in  the  algorithm 
and  it  produces  lower  misadjustment  than  the  Widrow-Hoff  LMS 
when  the  correlation  is  low  between  the  received  vector  and 
the  desired  response  (6:705). 


CHAPTER  3  :  ADAPTIVE  SYSTEM  MODELING 


3-1  INTRODUCTION 

This  chapter  will  outline  the  model  used  to  test 
and  simulate  the  Griffiths  adaptive  matched  filter.  First, 
an  overall  system  model  will  be  presented.  This  system 
model  will  then  be  expanded  and  broken  down  into  individual 
subsystems.  Each  subsystem  will  be  analyzed  separately. 

The  following  models  will  be  explained  in  this  chapter; 

1)  system  model 

2)  signal  and  composite  channel  filter  model 

3)  adaptive  matched  filter  model 

4)  Griffiths  algorithm  model 

5)  interference  and  noise  model 


3-2  A  SYSTEM  MODEL  FOR  ANALYSIS 

The  purpose  of  this  thesis  is  to  analyze  the 
performance  of  a  Griffiths  adaptive  matched  filter  over  a 
range  of  signal  inputs.  The  general  block  diagram  of  the 
system  model  is  shown  in  Figure  3-1.  The  input  signal  is 
8(k),  which  is  a  real  random  bit  stream  of  a  nonreturn-to- 
zero  level  (NRZ-L)  waveform  (11:79).  The  signal  8(k)  will 
then  pass  through  a  composite  channel  filter.  The  composite 
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interference 


noise 


Figure  3-1  General  Block  Diagram  of  the 
System  Model 


22 


channel  filter  will  represent  a  filter  that  both  pulse 

shapes  and  phase  distorts  the  original  bit  stream.  The  I 

composite  channel  pulse  shaping  filter  will  be  of  the  raised 
cosine  type  and  will  introduce  various  levels  of  phase 

distortion  (see  Appendix  B) .  The  composite  channel  filter  I 

finite  impulse  response  is  E(k) .  The  output  of  the 
composite  channel  filter  is  i(k) .  After  the  composite 

channel,  the  signal  will  have  noise  and  interference  added  I 

to  it.  The  signal  r(k) ,  which  has  channel  distortion  along 
with  noise  and  interference  added,  will  serve  as  input  to 
the  adaptive  matched  filter.  The  response  of  the  filter, 

Y(k),  to  various  inputs  (that  have  distortion,  noise,  and 
interference  added)  will  then  be  accomplished. 

The  signals  used  in  this  thesis  will  be  assumed  to 
be  narrowband.  This  implies  that  the  spectral  components  of 
the  signals  of  interest  are  confined  to  a  band  that  is  small 
compared  to  the  band  center  frequency,  or  carrier  frequency. 

In  general,  a  narrowband  representation  for  the 
input  waveform  can  be  expressed  as  (12:57) 

s(k)  =  ag(k)  cos[  p^k  +  es(k)  ]  (9) 

and  the  narrowband  representation  for  the  composite  channel 
filter  can  be  expressed  as 

p(k)  =  2  ap(k)  cos(  p^k  +  ep(k)  ]  (10) 
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where 


a(k)  =  amplitude  or  envelope  modulation 
0(k)  =  angle  modulation 
Pq  =  carrier  frequency 
k  =  the  discrete  time  index 

Note  that  in  equation  (10)  a  constant  factor  of  two  has  been 
added  for  convenience.  It  is  assumed  that  the  carrier 
frequency  is  much  greater  than  the  bandwidth  of  the  signal 
and  channel  filter. 

In  general  some  form  of  modulation  is  required  to 
transmit  information  across  a  communication  channel.  In 
this  thesis,  the  specific  form  of  modulation  will  not  be 
considered;  instead,  baseband  effects  will  be  analyzed  and 
the  complex,  low-pass,  baseband  representation  of  all 
seq\aences  will  be  used.  This  means  the  various  vectors  in 
Figure  3-1  can  be  defined  in  the  following  way 


E(k)  =  Ei(k)  +  j  EQ(k)  (11) 
i(k)  *  4i(k)  +  j  E^Ck)  (12) 
l(k)  =  4i(k)  +  j  4Q(k)  (13) 
n(k)  =  Bi(k)  +  j  BQ(k)  (14) 
r(k)  =  ri(k)  +  j  iQ(k)  (15) 

y(k)  =  Yi(k)  +  j  yQ(k)  (16) 
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where 


j  =  imaginary  operator 

subscript  I  =  the  in-phase,  or  I  channel,  component 
subscript  Q  =  the  quadrature,  or  Q  channel,  component 
k  =  the  iteration  or  clock  cycle 

Note  that  the  input  signal  £(k)  is  real  and  has  no 
imaginary  component. 


3-3  SIGNAL  AND  COMPOSITE  CHANNEL  FILTER  MODEL 

The  input  signal  is  fi{k),  which  is  a  real  random 
bit  stream.  For  instance,  s(k)  could  possibly  be  the 
information  content  of  a  radio  frequency  pulse  train 
generated  by  a  radar  or  a  digitally  coded  pulse  train  of  a 
radio  transmission.  s(k)  will  be  a  sampled  version  of  a 
NRZ-L  waveform.  Impulse  sampling  will  be  assumed  and  each 
sample  will  be  represented  by  a  Dirac  delta  function,  which 
is  also  known  as  an  impulse  function  (14:50).  The  bit 
stream  will  then  pass  through  a  composite  channel  filter 
which  will  pulse  shape  the  impulse  samples  and  introduce 
various  types  of  phase  distortion. 

The  composite  channel  filter  can  be  represented  as 
shown  in  Figure  3-2.  In  Figure  3-2  the  input  bit  stream  is 
passed  through  the  composite  channel  filter  and  the  output 
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Figure  3-2  Composite  Channel  Filter 
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is  two  signals,  Zj(k)  and  SQ(k) ,  which  represent  the  in- 
phase  and  quadrature  components  of  £(k) . 

The  transfer  function  of  the  composite  channel 
filter  can  be  expressed  as 

P(p)  =  A(3)  e‘3  *0)  (17) 

where 

PO)  =  the  channel  transfer  function 
A(p)  =  the  amplitude  response 
♦  O)  =  the  phase  response 

P  =  the  independent  frequency  variable 

In  order  to  generate  pulses  that  exhibit  some  form  of  phase 
distortion,  A(p)  is  set  equal  to  one  and  *(P)  is  defined 
the  following  way 

♦(P)  =  dQ  +  d^p  +  d2p^  +  +  d4P^  (18) 

where 

dQ  =  constant  phase  term 
d^  =  constant  delay  term 
d2  =  linear  delay  distortion  term 
d3  =  quadratic  delay  distortion  term 
d^  =  cubic  delay  distortion  term 


E-. 
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Note  that  more  terms  could  be  added  to  equation  (18)  to 
completely  describe  any  phase  characteristic.  "The  constant 
phase  term  and  the  constant  delay  term  alone  introduce  no 
distortion;  however,  they  play  an  important  part  in 
mathematically  describing  some  arbitrary  phase 
characteristic"  (17:9). 

If  the  composite  channel  filter  of  Figure  3-2  is 
implemented  as  a  raised  cosine  filter  with  roll-off  factor 
(see  Appendix  B)  equal  to  one  (see  Figure  3-3)  and  the 
spectrum  of  the  filter  is  given  by 


-a-  cos2 
2n 


0  <  1 p I  <  2n 


P(P)  =  < 


(19) 


0 


otherwise 


where 

ir/2n  =  a  scaling  factor  chosen  for  convenience 
2/T  =  the  signaling  rate 
T  =  ir/n 

2n  =  the  maximum  frequency  in  the  passband 
n  =  the  minimum  Nyquist  bandwidth 


then  the  components  of  the  composite  channel  filter  can  be 
expressed  as  (17:12) 

! 
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Figure^ 3-3  Raised  Cosine  Frequency  Spectrum 
with  Roll-Off  Factor  Equal  to  1 
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where 

=  normalized  time  variable  in  number  of  bit 
^  intervals 

Kq  =  constant  phase  component  expressed  in  radians 
=  constant  delay  component 
k2  =  linear  delay  component 
k3  =  quadratic  delay  component 
k^  =  cubic  delay  component 


Equations  (20)  and  (21)  give  the  composite  channel  filter 
impulse  response  for  the  in-phase  and  quadrature  components 
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The  outputs  of  the  composite  channel  filter  are  then  defined 
by  the  following  equations: 


ll(k)  =  a(k)  *  fii(k)  (22) 

Sg(k)  =  a(k)  *  EQ(k)  (23) 

By  selecting  different  values  of  kQ  thru  k^,  one  can 
generate  different  forms  of  phase  distortion  on  the  input 
bit  stream  s(k) . 

In  Figure  3-1,  the  system  model,  the  input  a 
needs  to  be  convolved  with  the  channel  response.  To 
facilitate  the  convolution  of  these  functions  the  fast 
Fourier  transform  (FFT)  implementation  of  convolution  will 
be  used.  Figure  3-4  shows  how  the  FFTs  will  be  used  to 
perform  the  linear  convolution.  The  input  sequences  need  to 
be  augmented  with  zeros  so  that  the  required  FFT 
multiplication  results  in  a  linear  convolution  (7:270). 


3-4  ADAPTIVE  MATCHED  FILTER  MODEL 

The  adaptive  matched  filter  will  be  implemented  as 
a  tapped-delay-line.  Figure  2-4  shows  the  implementation. 

In  general,  there  are  a  total  of  N  samples  per  symbol 
waveform,  where  N  is  an  integer.  In  this  thesis  N=8  and 
the  vector  length  of  C(k)  will  be  8.  Y(k)  is  sampled  at 
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Figure  3-4  FFT  Implementation 


time  t=T,  which  is  also  the  time  required  to  justify  £ 
within  the  matched  filter.  Thus  one  sample  of  Y(k)  occurs 
for  each  vector  input  £. 

The  weight  vector  if(k) ,  which  is  composed  of  Wj^ 
thru  Wjj,  is  adaptive  so  that  the  individual  weights  can  be 
updated  through  the  adaptive  algorithm. 


3-5  GRIFFITHS  ALGORITHM  MODEL 

The  complex  Griffiths  algorithm  was  derived 
previously  and  is  given  by 


w(k+l)  =  w(k)  +  2n  £  -  2m.  Y(k)  r*(k) 
where  5  =  B{  D(k)  £*(k)  > 


(24) 


All  terms  in  the  algorithm  are  known  apriori  or 
can  be  measured  directly  from  observed  data.  The  P-vector 
requires  further  elaboration  to  gain  a  full  understanding. 
The  complex  P-vector  is  defined  to  be 

£  =  E{  D(k)  £*(k)  )  (25) 


where 


E  =  the  statistical  expectation  operator 
D(k)  =  the  complex  desired  signal 
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X*(k)  =  the  complex  conjugate  of  the  received  vector 

D(k) ,  for  a  communication  application,  is  a  complex  constant 
determined  by  the  initial  receiver  design  and  signal 
constellation.  With  this  assumption  the  P-vector  can  be 
written  as 

£  =  D(k)  B{  l*(k)  )  (26) 

Assuming  now  that  we  have  no  interference  and  using  notation 
consistent  with  Figure  3-1,  the  P-vector  can  be  written  as 

£  =  D(k)  E{  ii(k)  -  j  r^Ck)  }  (27) 

£  =  D  E{  (ij  +  Dj)  -  j  (Sq  +  Sq)  )  (28) 

where  the  discrete  time  index,  k,  has  been  dropped  for 
convenience.  Assuming  that  the  noise  is  a  zero  mean 
process,  the  P-vector  becomes 

£  =  D  [  E{  ij  )  -  j  E{  Sg  )  ]  (29) 

£  =  (Dj  +  j  Dq)  [  E{  Zj  )  -  j  E(  ^  )  ]  (30) 

The  values  Ej  and  Zg  are  outputs  of  the  composite 
channel  filter.  It  can  be  shown  (13:172)  that  for  a  causal, 
time- invariant  linear  filter,  the  mean  of  the  output  is 
given  by 
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Diy(k)  = 


(31) 


where 

=  the  mean  of  the  output 
mjj(k)  =  the  mean  of  the  input 
h(k)  =  a  channel's  finite  impulse  response 


Therefore,  the  expectations  given  in  equation  (30)  can  be 
written  as 


B{2l) 


(32) 


(33) 


Equations  (24),  (30),  (32),  and  (33)  define  the 
Griffiths  algorithm  as  implemented  in  the  simulation. 


3-6  INTERFERENCE  AND  NOISE  MODEL 

Both  interference  and  noise  will  be  modeled  as 
complex,  low-pass,  baseband  signals.  Therefore,  both  can  be 
defined  as 


l(k)  =  ij.(k)  +  j  ig(k)  (34) 

B(k)  =  Bj;(k)  +  j  ngCk)  (35) 
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The  noise  process  will  be  approximated  as  an 
additive  white  gaussian  process  with  zero  mean.  In  order  to 
calculate  one  element  of  the  noise  matrix  the  following 
equation  will  be  used 


n(k) 


(36) 


where 

n(k)  =  a  gaussian  noise  value 

o  =  the  standard  deviation  of  the  noise  process 
Uj^  =  a  uniformly  distributed  random  variable 
between  0  and  1 

Equation  (36)  illustrates  how  the  noise  is  generated  in  the 
computer  simulation.  The  basis  for  equation  (36)  lies  in 
the  central  limit  theorem.  This  theorem  states  that  the  sum 
of  N  independent,  identically  distributed  random  variables 
becomes  a  Gaussian  distribution  as  N  becomes  large 
(18:243).  In  ecjuation  (36)  N=12,  and  the  Gaussian 

approximation  is  based  on  the  sum  of  12  uniform  random 
variables.  If  the  noise  is  also  assumed  to  be  ergodic  (see 
Appendix  B)  then  the  standard  deviation  is  related  to  the 
noise  power  by 

o  =’\Jnoise  power  (37) 
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Equation  (37)  enables  one  to  specify  a  noise  power  and 
generate  noise  at  this  power  level. 

The  interference  will  also  be  additive.  The 
interference  will  be  narrowband  and  be  in  the  same  frequency 
band  as  the  signal  of  interest;  therefore,  i(k)  takes  the 
form  of  equation  (34)  .  The  in-phase  and  quadrature 


components  can 

be  represented  as 

il(k)  =• 

^2  JPWR  cos[  2if  Fj  Tj  k  + 

2it  cl] 

(38) 

iQ(k)  =• 

^2  JPWR  sin[  2it  Fj  Tj  k  + 

2it  cl] 

(39) 

where 

JPWR 

=  the  jamming  power 

Fj 

=  the  frequency  of  the  jammer 

Tj 

=  the  period 

k 

=  the  discrete  time  index 

cl 

=  the  phase  component 

Equations  (38)  and  (39)  define  the  interference  model  as 
used  in  the  simulation. 


37 


J 


3-7  GENERAL  SIMULATION  MODEL 

The  purpose  of  this  thesis  is  to  analyze  the 
performance  of  a  complex  Griffiths  adaptive  matched  filter 
for  various  input  signals.  The  software  chosen  for  this 
analysis  is  MathCAD  developed  by  Mathsoft,  Inc.  A  complete 
listing  of  the  simulation  software  program  is  found  in 
Appendix  A.  The  general  procedure  for  the  simulation  is  as 
follows 


STEP  1: 
STEP  2: 
STEP  3: 
STEP  4: 
STEP  5: 

STEP  6; 

STEP  7; 


define  s (k) . 

define  Ei(k)  and  Eq(I^)  • 
define  D(k) ,  the  desired  response, 
define  ijCk),  igCk),  ajCk),  Bq(Ic)  . 
input  the  defined  signals  and  run  the 
simulation. 

compare  the  output  Y(k)  against  the 
various  defined  input  signals 
compare  the  performance  of  the  adaptive 
filter  against  the  various  defined  input 
signals. 


Note  that  the  simulation  requires  the  user  to  define  the 
various  input  sequences. 


38 


Throughout  the  development  of  the  thesis  model 
various  assumptions  were  made  concerning  the  filter  and 
model.  All  assumptions  are  repeated  here  for  compactness: 

1)  input  signal  is  at  least  cyclostationary  and 
sampled  at  or  above  the  Nyquist  rate.  The  input  signal  is 
real,  has  odd  length,  and  is  antipodal. 

2)  synchronization  is  maintained  during  all 
operations . 

3)  all  signals  are  modeled  as  narrowband,  complex, 
low-pass,  baseband  signals. 

4)  only  filtering  operations  are  considered. 

5)  all  filters  are  stable,  causal,  FIR  filters. 

6)  channel  filter  is  linear,  time-invariant. 

7)  the  adaptive  filter  is  implemented  as  a  FIR 
structure  and  uses  closed  loop  adaptation. 

8)  complex  Griffiths  algorithm  is  implemented 

9)  noise  and  interference  are  both  additive. 
Interference  has  a  zero  mean. 

10}  noise  is  approximated  as  a  zero  mean, 
independent,  additive  white  gaussian  noise  process. 

11)  noise,  interference,  and  input  signal  are 


ergodic. 


CHAPTER  4  :  ANALYSIS  OF  RESULTS 


4-1  INTRODUCTION 

This  chapter  analyzes  performance  of  the  Griffiths 
adaptive  matched  filter  over  a  range  of  signal  inputs.  All 
notation  will  be  consistent  with  that  of  Figures  3-1  and  3- 
2.  First,  several  performance  criteria  are  established. 

This  discussion  will  lay  the  fundamental  framework  for  the 
jamming  analysis  that  follows.  The  jamming  analysis  will  be 
conducted  over  a  range  of  signal  distortion  levels  and 
jamming  to  signal  power  ratios  (JSR) . 


4-2  PERFORMANCE  CRITERIA 

4-2-1  POWER  CRITERIA 

Throughout  the  analysis  the  terms  signal  power, 
noise  power,  and  jamming  power  will  be  used.  Signal  power 
will  be  the  power  in  A(k) ,  noise  power  will  be  the  power  in 
a(k) ,  and  jamming  power  will  be  the  power  in  i(k) . 


Throughout  the  analysis  the  following  constants  will  be 
assumed 

U  =  •  01 

D  =  1  +  j 

signal  power  =  .375 
noise  power  =  0 

s(k)  =  the  same  sequence  of  binary 
data 

The  reason  for  setting  the  noise  power  to  zero  is  to 
explicitly  analyze  the  effects  of  jamming  and  delay 
distortion  on  the  filter. 

The  convergence  rate  constant,  n,  is  set  equal  to 
.01  so  that  jamming  power  levels  can  be  set  equal  to  -10  DB, 
0  DB,  and  10  DB.  In  general,  the  convergence  rate  constant 
should  satisfy  (8:103) 

1 

U  <  -  (40) 

(number  of  tap  weights) (received  power) 

rearranging  and  substituting 

1 

received  power  <  —  ■  ■■  .■  —  (41) 

(8)  (.01) 


received  power  <  12.5 


(42) 


since  the  transmitted  signal  power  is  .375,  the  combined 
noise  and  jamming  power  must  be  no  greater  than  12.125. 

This  implies  that  the  maximum  JSR  is  approximately  15  DB. 
With  the  convergence  rate  constant  equal  to  .01,  analysis  of 
JSRs  equal  to  -10  DB,  0  DB,  and  10  DB  is  possible. 

s(k)  will  be  a  real  bit  stream  of  length  31. 
s(k)  will  be  generated  as  a  pseudonoise  sequence  of  length 
31  and  take  the  following  form  in  binary  data 

1000  0100  1011  0011  1110  0011  0111  010 

The  sequence  for  s(k)  was  generated  from  a  5-stage  linear 
feedback  shift  register  as  shown  in  Figure  4-1  with  an 
initial  shift  register  state  of  00001.  The  secjuence  for 
s(k)  is  of  maximal  length  and  can  be  shown  to  obey  the 
properties  of  pseudonoise  sequences  such  as  balance,  run, 
and  correlation  (11:546).  A  plot  of  the  normalized 
autocorrelation  function  is  shown  in  Figure  4-2. 

4-2-3  AVERAGE  ERROR  CRITERION 

For  this  thesis  the  error  will  be  defined  as 

e(k)  =V[  Re{D)  -  Re(Y)  ]2  +  [  lm(D)  -  Im(Y)  ]2(43) 

where 

e(k)  =  the  error  signal 
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Figure  4- 


1  5-Stage  Linear  Feedback 
Shift  Register 


I 


I 
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V# 


D  =  the  desired  signal 
Y  =  the  actual  filter  output 
PNjjjt  ~  transmitted,  +1  or  -1 

e(k)  is  simply  the  distance  between  D  and  Y  If 

e(k)  is  the  actual  error  at  a  specific  adaptation  cycle, 
then  the  average  error  will  be  defined  as 

1  185 

r  =  n  E  e(Jt)  (44) 

N  k=86 


where 

N  =  the  number  of  adaptation  cycles  (100) 

r  *  the  average  error 

The  average  error  will  represent  the  average  of  the  error 
over  the  last  100  adaptations.  This  is  done  to  minimize 
transient  effects. 

To  make  a  comparison  to  other  input  signals  the 
number  of  adaptations  will  be  fixed  at  186.  186  was  chosen 

because  after  this  many  adaptations  the  filter  has  settled 
down  well  into  its  steady  state  for  ii=.01  (for  instance,  see 
Figure  4-7) . 

4-2-4  DEFINING  Dj^^k)  AND 

The  in-phase  and  quadrature  response  of  the 
composite  channel  filter  are  defined  by  equations  (20)  and 


45 


(21) .  For  this  thesis  we  will  be  concerned  with  four  cases. 
These  cases  are 


CASE  I  :  no  distortion 
CASE  II  :  linear  delay  distortion 
CASE  III  :  quadratic  delay  distortion 
CASE  IV  ;  cubic  delay  distortion 

The  representations  for  Ei(Jt)  and  EgCk)  are  given  for 
each  of  the  four  cases  in  Figures  4-3  thru  4-6.  The  amount 
of  delay  distortion  is  set  at  k/T  =  2. 


4-2-5  DEFINING  THE  P-VECTOR 

The  P-vector  was  defined  in  equation  (25) .  One  of 
the  requirements  of  the  Griffiths  algorithm  is  that  the  P- 
vector  is  determined  apriori  (4:1699).  For  this  thesis  the 
P-vector  will  be  determined  from  Figure  4-3  for  the  case  of 
no  distortion.  With  this  assumption  the  P-vector  becomes 


P  =  D 


0 

.17 

.5 

.849 

1 

.849 

.5 

.17 


Note  that  the  P-vector  is  determined  from  the  main  lobe  of 
the  in-phase  channel  and  consists  of  8  sample  values.  The 
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Figure  4-4  Case  II.  Linear  Delay  Distortion. 
Calculated  from  Equations 
(20)  and  (21)  with  T=.5,  kQ=0, 

}C2^“2 /  ^2* 1/  0/  k^~0. 
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Figure  4-5  Case  III.  Quadratic  Delay  Distortion. 
Calculated  from  Equations 
(20)  and  (21)  with  T=.5,  )Cq=0, 

)C2^=2 ,  K2~0,  k^=0 • 
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Figure  4- 


6  Case  IV.  Cubic  Delay  Distortion. 
Calculated  from  Equations 
(20)  and  (21)  with  T=.5, 

^2^—2,  f  )C2“0^  )c^“l. 
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P-vector  has  length  8  because  this  is  the  length  of  the 
filter. 


4-3  JAMMING  RESULTS 

This  section  will  present  the  results  of  the 
simulation.  The  first  section  will  present  the  effects  of 
jamming  power  and  jamming  freguency  on  the  filter  output. 

The  second  section  will  present  the  effects  of  jamming  phase 
on  the  filter  output.  Subsequent  sections  will  compare  CW 
jamming  and  noise  jamming  and  make  comparisons  to  the  LMS 
algorithm. 


4-3-1  EFFECTS  OF  JAMMING  POWER  AND  JAMMING 
FREQUENCY  (PHASE  CONSTANT) 

The  first  set  of  results  is  that  for  zero  jamming 
power  and  zero  noise  power.  These  results  are  presented  in 
Figures  4-7  thru  4-11.  The  actual  values  for  the  average 
error  are  presented  in  Tables  1  thru  3  for  a  JSR  =  Note 

that  pulses  that  experience  linear  delay  distortion  will  not 
enable  the  algorithm  to  converge.  This  will  be  explained 
further  in  the  COMPARISONS  section. 

The  results  of  the  simulation  dealing  with  jamming 
power  and  jamming  frequency  (phase  constant)  are  presented 
in  Tables  1  thru  3.  A  graphical  presentation  of  these 
results  are  presented  in  Figures  4-12  thru  4-14. 
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NO  JAMMING 

No  Distortion 


Figure  4-7  No  Jamming/No  Distortion 


NO  JAMMING 

Linear  Delay 
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Jamming/Linear  Delay  Distortion 


NO  JAMMING 

Quadratic  Delay 
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Figure  4-9  No  Jamming/Quadrat ic  Delay  Distortion 


NO  JAMMING 

Cubic  Delay 


Figure  4-10  No  Jamming/Cubic  Delay  Distortion 


NO  JAMMING 


HOHHa  aovnaAV 
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No  Distortion  Quadratic  Delay  Cubic  Delay 


TABLE  1  Average  error  of  filter  output  for 
an  input  that  has  no  distortion 


JSR  (DB) 

JAMMING 
CYCLES  per 

PN  LENGTH 

AVERAGE 

ERROR 

(r) 

—00 

- 

.002 

-10 

1 

.332 

0 

1 

.371 

10 

1 

.744 

-10 

2 

.26 

0 

2 

.296 

10 

2 

.751 

-10 

8 

.2 

0 

8 

.237 

10 

8 

.896 

-10 

31 

.  079 

0 

31 

.  127 

10 

31 

.799 

Note;  Noise  power  is 

zero  and  the  phase 

of  the  jammer  is 

zero 
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TABLE  2 

Average  error  of  filter  output  for 
an  input  that  has  quadratic  delay 
distortion 

1 

i 

JSR  (DB) 

JAMMING 

AVERAGE 

CYCLES  per 

ERROR 

- 

PN  LENGTH 

(r) 

1 

-00 

- 

.038 

-10 

1 

.319 

i 

0 

1 

.416 

li 

10 

1 

.623 

-10 

2 

.28 

0 

2 

.  376 

1 

10 

2 

.647 

-10 

8 

.318 

J 

0 

8 

.425 

1 

10 

8 

.911 

Note;  Noise  power 

is 

zero  and  the  phase  of 

the 

jammer  is 

J 

zero 
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TABLE  3  Average  error  of  filter  output  for 

an  input  that  has  cubic  delay  distortion 


JSR  (DB) 

JAMMING 
CYCLES  per 

PN  LENGTH 

AVERAGE 

ERROR 

(r) 

-00 

- 

.083 

-10 

1 

.369 

0 

1 

.403 

10 

1 

.76 

-10 

2 

.  304 

0 

2 

.352 

10 

2 

.751 

-10 

8 

.262 

0 

8 

.325 

10 

8 

.888 

Note:  Noise  power  is 

zero  and  the  phase 

of  the  jammer  is 

zero 


10  DB 


CYCLES  per  PN  LENGTH 

No  Distortion  Quad  Delay  Cubic  Delay 


JSR  =  0  DB 


CYCLES  per  PN  LENGTH 

No  Distortion  Quad  Delay  Cubic  Delay 

Figure  4-13  o  DB  Plot  of  Jamming  Results  (Phase=0) 


10  DB 


} 

HOHHa  aovHaAV 
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CYCLES  per  PN  LENGTH 

No  Distortion  Quad  Delay  Cubic  Delay 

‘igure  4-14  10  DB  Plot  of  Jamming  Results  (Phase=0) 


The  effects  of  jamming  phase  at  0  OB  and  -10  DB 
are  presented  in  Tables  4  thru  9.  Figures  4-15  thru  4-20 
are  a  graphical  presentation  of  these  results. 

4-3-3  EFFECTS  OF  NOISE  ON  FILTER 

In  order  to  compare  between  noise  jamming  and 
continuous  wave  jamming,  a  one-time  run  was  made  on  a  signal 
that  has  no  distortion  at  various  signal-to-noise  ratios 
(SNR) .  These  results  are  presented  in  Table  10.  Actual 
error  output  due  to  various  SNRs  is  graphically  depicted  in 
Figures  4-21  thru  4-23. 

4-3-4  LMS  ALGORITHM 

The  purpose  of  this  thesis  has  been  to  analyze  an 
adaptive  matched  filter  that  uses  the  Griffiths  algorithm. 
For  comparison  several  selected  runs  were  repeated  using  the 
LMS  algorithm.  The  LMS  algorithm  is  implemented  as 

w(k+l)  =  w(k)  +  2nPNj,.jD  r*(k)  -  2n  Y(k)  r*(k)  (45) 


where  all  terms  have  been  defined  previously.  The  LMS  runs 
are  presented  in  Figures  4-24  thru  4-27  and  Table  11. 


TABLE  4  Average  error  of  filter  output  for 
different  phase  constants  (cl) . 

JSR=0  DB  and  Frequency  =  1  cycle  per 
PH  Length. 


No  Distortion  cl  r 


Quadratic  Delay 

0 

.25 

.5 

.75 

.371 

.263 

.273 

.386 

0 

.416 

.25 

.342 

.5 

.313 

.75 

.393 

Cubic  Delay 

0 

.403 

.25 

.349 

.5 

.434 

.75 

.513 

f 


TABLE  5  Average  error  of  filter  output  for 
different  phase  constants  (cl) . 

JSR=0  DB  and  Frequency  =  2  cycles  per 
PN  Length. 


No  Distortion  cl  r 


Quadratic  Delay 

0 

.25 

.5 

.75 

.296 

.269 

.367 

.4 

0 

.  376 

.25 

.334 

.5 

.363 

.75 

.431 

Cubic  Delay 

0 

.352 

.25 

.417 

.5 

.504 

.75 

.441 
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TABLE  6  Average  error  of  filter  output  for 
different  phase  constants  (cl) . 

JSR=0  DB  and  Frequency  =  8  cycles  per 
PN  Length. 


No  Distortion 

cl 

r 

0 

.237 

.25 

.309 

.5 

.387 

.75 

.325 

Quadratic  Delay 

0 

.425 

.25 

.326 

.5 

.378 

.75 

.529 

Cubic  Delay 

0 

.325 

.25 

.44 

.5 

.455 

.75 

.318 
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TABLE  7  Average  error  of  filter  output  for 
different  phase  constants  (cl) . 

JSR=-10  DB  and  Frequency  =  1  cycle  per 
PN  Length. 


No  Distortion  cl  r 


0  .332 
.25  .241 
.5  .257 
.75  .341 


Quadratic  Delay 


0 

25 

5 

75 

.319 

.255 

.251 

.315 

0 

.369 

25 

.291 

5 

.313 

75 

.  385 

Cubic  Delay 


L 


TABLE  8  Average  error  of  filter  output  for 

different  phase  constants  (cl) .  _ 
JSR=~10  DB  and  Frequency  =  2  cycles  per  I 
PN  Length. 


No  Distortion  cl  F 


Quadratic  Delay 

0 

.25 

.5 

.75 

.26 

.234 

.327 

.344 

0 

.  28 

.25 

.241 

.5 

.295 

.75 

.333 

Cubic  Delay 

0 

.  304 

.25 

.3 

.5 

.371 

.75 

.382 

H 


I 
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TABLE  9  Average  error  of  filter  output  for 
different  phase  constants  (cl) . 

JSR=-10  DB  and  Frequency  =  8  cycles  per 
PN  Length. 


No  Distortion  cl  r 


Quadratic  Delay 

0 

.25 

.5 

.75 

.2 

.275 
.  34 
.279 

0 

.318 

.25 

.263 

.5 

.344 

.75 

.  386 

Cubic  Delay 

0 

.262 

.25 

.  341 

.5 

.  387 

.75 

.311 
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PHASE  ERROR  AT 

1  CYCLE  per  PN  LENGTH 


®  m  Tf  P7  w  ^  o 

d  d  d  d  d  d  d 

JOJJ3  aSBJSAV 
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Phase  Constant  cl 

No  Distortion  -*■  Quad  Delay  Cubic  Delay 

Figure  4-15  Average  Error  Due  to  Phase  (0  DB  and 
1  Cycle  per  PN  Length) 


PHASE  ERROR  AT  0  DB 

2  CYCLES  per  PN  LENGTH 


J0JJ3  99iej3Ay 
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Phase  Constant  cl 

Distortion  Quad  Delay  Cubic  Delay 

jure  4-16  Average  Error  Due  to  Phase  (0  DB  and 
2  Cycles  per  PN  Length) 


PHASE  ERROR  AT  0  DB 

0  CYCLES  per  PN  LENGTH 


Phase  Constant  cl 

Distortion  Quad  Delay  Cubic  Delay 

jure  4-17  Average  Error  Due  to  Phase  (0  DB  and 
8  Cycles  per  PN  Length) 


PHASE  ERROR  AT  -10  DB 

1  CYCLE  per  PN  LENGTH 
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Phase  Constant  cl 

Distortion  Quad  Delay  Cubic  Delay 

ire  4-18  Average  Error  Due  to  Phase  (-10  DB 
1  Cycle  per  PN  Length) 


PHASE  ERROR  AT  -10  DB 

2  CYCLES  per  PN  LENGTH 
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Phase  Constant  cl 

No  Distortion  Quad  Delay  Cubic  Delay 

Figure  4-19  Average  Error  Due  to  Phase  (-10  DB  and 
2  Cycles  per  PN  Length) 


PHASE  ERROR  AT  -10  DB 

8  CYCLES  per  PN  LENGTH 
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Phase  Constant  cl 

No  Distortion  Quad  Delay  Cubic  Delay 

Figure  4-20  Average  Error  Due  to  Phase  (-10  DB  and 
8  Cycles  per  PN  Length) 


TABLE  10  Average  error  of  filter  output 

for  various  SNRs  (No  Distortion) 


No  Distortion 

SNR  =  10  DB 


HOHHa  ivniov 
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ADAPTATIONS 

Figure  4-21  Actual  Error  Due  to  Noise  (SNR  =  10  DB) 


-22  Actual  Error  Due  to  Noise 


No  Distortion 


ADAPTATIONS 


LMS  ALGORITHM 

Linear  Delay  Distortion  (No  Jamming) 


ADAPTATIONS 

Figure  4-24  Linear  Delay  Distortion,  No  Jamming, 


LMS  ALGORITHM 

lie  Delay  Distortion  (No  Jamming) 


ADAPTATIONS 

Figure  4-25  Cubic  Delay  Distort 
LMS  Algorithm 


LMS 


CD  cvi  o  cq 

^  t-H  tH  O 


CD 

o' 


HOHHa  ivniov 
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0.4' 


LMS  ALGORITHM  (No  Distortion) 

JSR=10  DB,  0  CYCLES  per  PN  LENGTH 


OCOCD’^CMOCqtp'^WO 

HOHHa  ivniDV 
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ADAPTATIONS 

Figure  4-27  LMS  Algorithm  (No  Distortion) ,  JSR  =  10  DB 
8  Cycles  per  PN  Length 


TABLE  11  Comparisons  of  Griffiths  and  LMS 
algorithms  for  selected  runs 


r 


UlS 

GRIFFITHS 

073 

DIVERGES 

017 

.083 

203 

.371 

RUN 


Linear  Delay 

(No  Noise,  No  Jamming) 

Cubic  Delay 

(No  Noise,  No  Jamming) 

No  Distortion,  JSR=0  DB 
1  Cycle  per  PN  Length 


124 


,896 


No  Distortion,  JSR=10  DB 
8  Cycles  per  PN  Length 


4-4  CQWPARISQNg 

This  section  will  compare  various  results  obtained 
in  Section  4-3. 

4-4-1  COMPARISONS  OF  NO  NOISE  AND  NO  JAMMING 
RESULTS 

One  of  the  most  significant  results  is  that  the 
Griffiths  algorithm  diverges  when  pulses  that  experience 
linear  delay  distortion  are  processed  by  the  filter.  This 
divergence  occurs  even  though  the  received  signal  power  is 
well  within  the  limits  established  by  eguation  (42) . 

The  linear  delay  distortion  run  of  Figure  4-8  was 
executed  again  for  (1)  =  .001,  adaptations  =  930  and  (2) 

«  .01,  adaptations  =  1085.  The  results  were  similar  to 
Figure  4-8  again  —  the  algorithm  was  not  converging. 

By  looking  at  Table  12  it  is  shown  that  for  the 
case  of  linear  delay  distortion  the  power  in  the  main  lobe 
(i.e.  that  portion  the  matched  filter  uses  to  integrate 
upon)  is  the  least  compared  to  the  other  types  of 
distortions.  This  means  linear  delay  distortion  causes  the 
most  spreading  of  the  pulse.  Also,  linear  delay  distortion 
is  the  least  correlated  with  the  P-vector.  The  complex 
correlation  coefficient  is  defined  as  (21:291) 

Cov(X, Y) 

B  =  -  (46) 
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TABLE  12  Comparisons  of  various  distortion 
levels  (no  jamming  or  noise) 


Distortion  Level 

Correlation 

B 

Coefficient 

|b| 

%  of  power 
in  main 
lobe 

No  Distortion 

1 

1 

100% 

Linear  Delay 

.773  -  .345j 

.847 

72.8% 

Quad  Delay 

.959 

.959 

98.4% 

Cubic  Delay 

.99  ~  .05j 

.991 

98.1% 
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where 


B  =  the  complex  correlation  coefficient 
X,  Y  =  complex  vectors 
Cov  =  the  covariance 

o  =  the  standard  deviation  of  one  complex  vector 

Table  12  shows  that  linear  delay  distortion  is  the  least 
correlated  with  the  P-vector  and  has  the  most  pulse 
spreading.  These  results  show  that  a  combination  of  (1) 
correlation  with  the  P-vector  and  (2)  pulse  spreading  have  a 
direct  impact  on  the  convergence  of  the  algorithm  even 
though  received  power  is  within  limits. 

The  LMS  algorithm  is  able  to  handle  linear  delay 
distortion.  This  is  shown  in  Figure  4-24.  Also,  the  LMS 
algorithm  produces  lower  average  error  than  the  Griffiths 
algorithm  for  cubic  delay  distortion  (r=.083  for  Griffiths 
and  r=.017  for  LMS), 

From  the  results  of  no  noise/no  jamming  it  can  be 
seen  that  the  LMS  adaptive  matched  filter  is  superior  to  the 
Griffiths  adaptive  matched  filter  in  handling  delay 
distortion  with  PN  sequences.  The  LMS  is  superior  because 
it  uses  a  stored  reference  which  represents  the  desired 
signal  exactly. 


NG  RESULTS 


The  results  of  Tables  1  thru  3  show  that  by 
increasing  jamming  power  for  a  given  frequency  the  average 
error  will  increase. 


4-4-2B  Effects  of  Jamming  Frequency /Phase 
The  Griffiths  filter  responds  well  to  frequency 
and  phase  changes.  The  average  error,  in  general,  changes 
with  frequency  and  phase  changes.  The  actual  magnitude  and 
direction  of  change  can  be  found  from  Tables  1  thru  9.  The 
Griffiths  filter  is  able  to  converge  to  a  steady  state  with 
frequency  and  phase  changes  in  a  jammer. 


4-4-2C  Effects  of  Noise  Jamming 
Increasing  the  power  of  the  noise  will  result  in  a 
larger  average  error.  Noise  jamming  produces  a  higher 
average  error  for  JSRs  of  10  DB  and  0  DB  compared  to  a  CW 
jammer.  At  -10  DB  JSR,  on  the  average,  a  CW  jammer  produces 
a  higher  average  error. 


4-4-2 D  COMPARISONS  TO  LMS 
The  selected  LMS  runs  are  compared  to  the  same 
Griffiths  runs  in  Table  11.  Based  on  these  results  the  LMS 
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adaptive  matched  filter  is  superior  to  the  Griffiths 
adaptive  matched  filter  for  PN  sequences  and  CW  jamming. 
The  LMS  is  superior  because  it  has  a  desired  signal 
available  in  the  form  of  a  stored  reference. 


CHAPTER  5  t  CONCLUSION 


5-1  SUMMARY 

This  thesis  presented  a  CW  and  noise  jamming 
analysis  of  an  adaptive  matched  filter  that  (1)  has  a  PN 
sequence  as  an  input  and  (2)  uses  the  Griffiths  algorithm. 
The  Griffiths  filter  is  shown  to  converge  over  JSRs  of  -10 
DB,  0  DB,  and  10  DB  and  frequencies  of  1,  2,  and  8  cycles 
per  PN  length.  Also,  the  filter  is  shown  to  converge  over 
jamming  phases  of  0,  .Btt,  it,  and  1.5ir  radians. 

The  Griffiths  filter  is  shown  to  converge  for 
raised  cosine  pulses  that  experience  quadratic  delay 
distortion  and  cubic  delay  distortion.  For  pulses  that 
experience  linear  delay  distortion,  the  filter  is  shown  to 
diverge  even  though  the  convergence  rate  constant  is  well 
within  limits.  This  divergence  is  due  to  the  lack  of 
correlation  between  the  input  signal  and  the  P-vector.  This 
divergence  for  pulses  that  experience  linear  delay 
distortion  is  a  disadvantage  of  using  the  Griffiths  adaptive 
matched  filter.  Throughout  the  analysis  the  P-vector  is 
determined  apriori  and  held  constant. 

The  Griffiths  adaptive  matched  filter  is  shown  to 
converge  for  both  CW  jamming  and  noise  jamming.  Noise 
jamming  is  more  effective  than  CW  jamming  in  the  0  DB  and  10 
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DB  JSR  ranges.  CW  jamming,  on  the  average,  is  more 
effective  in  the  -10  DB  JSR  range. 

A  comparison  is  made  between  the  Griffiths 
adaptive  matched  filter  and  an  adaptive  matched  filter  that 
uses  the  LMS  algorithm.  The  LMS  adaptive  matched  filter  is 
shown  to  be  superior  to  the  Griffiths  adaptive  matched 
filter  when  a  stored  reference  of  the  desired  signal  is 
available  for  use  with  the  LMS. 


5-2  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  following  are  areas  recommended  for  further 

study. 

1)  The  P-vector  was  held  constant  during  the 
analysis.  Further  research  should  be  conducted  in  the  area 
of  updating  the  P-vector  with  each  new  processing  block 

(i.e.  determine  the  P-vector  from  received  data)  to  see  if  ! 

performance  can  be  improved. 

2)  The  effects  of  high  jamming  frequencies  need  to 

be  explored  on  the  Griffiths  adaptive  matched  filter.  I 

3)  The  desired  signal  was  a  constant  in  this 

thesis.  Further  research  should  be  conducted  with  a  non¬ 
constant  desired  signal  to  see  if  performance  can  be  I 

improved . 

4)  It  was  shown  that  pulses  that  experience  linear 

delay  distortion  cause  the  filter  to  diverge  even  though  the  j 
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convergent  rate  constant,  \i,  is  within  limits.  This 
phenomenon  should  be  studied  further.  A  more  detailed 
correlation  analysis  concerning  the  received  signal  with  the 
P-vector  should  be  accomplished  in  order  to  explain  this 
divergence  more  thoroughly. 


I 


1 


i 
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Appendix  A  : 


Software  Simulation  Programs 


This  appendix  lists  the  computer  programs  used  in  the 
simulation.  All  programs  are  MathCAD  programs.  User 
instructions  for  MathCAD  can  be  found  in  reference  [16]. 
This  appendix  contains  the  following  programs: 

1)  JAM.MCD 

2)  PNCODE.MCD 

3)  REVEC.MCD 

4)  PVEC.MCD 

5)  GRIFFPV.MCD 

6)  CONTPV.MCD 

7)  NOISE. MCD 
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kmax  -  Number  of  rows  in  jamming  matrix.  Corresponds  to 
number  of  samples  per  symbol. 

IMAX  -  Number  of  columns  in  jamming  matrix.  Corresponds 
to  the  number  of  adaptation  cycles. 

VMAX  -  Number  of  elements  in  jamming  matrix. 

k,  ITER,  VLEN  -  Counting  variables. 

JPWR  -  Jamming  power. 

FJ  -  Jamming  frequency. 

T  -  Jamming  period. 

cl  -  Jamming  phase. 

il  -  In-phase  vector. 

iQ  ”  Quadrature  vector. 

JAMI  -  In-phase  matrix. 

JAMQ  -  Quadrature  matrix. 

PWRI  -  Power  in  the  in-phase  components. 

PWRQ  -  Power  in  the  quadrature  components. 

FILES  USED; 

JAMI.dat  -  data  file  that  stores  the  in-phase  values. 
JAMQ.dat  -  data  file  that  stores  the  quadrature  values. 
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^3t> 


■  ' — 

0 

JAMI 

:=  il 

VLEN 

k, ITER  k+ (kmax+l)  ■  ITER 

JAMQ 

:=  iQ 

k ,  ITER  k+  ( kinax+ 1 )  ■  ITER 

WRITEPRN 

■jAMI 

dat 

:=  JAMI 

WRITEPRN 

JAMQ 

dat 

:=  JAMQ 

>WRI  :=  — 

1  — 1  2 

-  >  lil  1 

VMAX  +  1  ^ 

— 1  1  VLENj 

VMAX  +  1 


NAME  OF  PROGRAM;  PNCODE.MCD 

PURPOSE ;  Generates  values  for  2j(k)  and  ZgCk) . 

SELECTED  VARIABLE  IDENTIFICATION; 

VMAX  -  FFT  vector  length. 

BMAX,  kmax  -  Number  of  rows.  Corresponds  to  number  of 
samples  per  symbol. 

TMAX,  IMAX  -  Number  of  columns.  Corresponds  to  number 
of  adaptations. 

PN  -  Pseudonoise  sequence. 

BIT  -  Column  vector  which  gives  length  and  definition 
of  the  bit. 

T  -  Column  vector  used  in  FFT  calculations. 

BITS  -  Matrix  of  pseudonoise  bits. 

PULSI,  PULSQ  -  I  and  Q  channels  of  the  pulse  shaping 
impulse  response. 

PI  thru  P5  -  FFT  processing  arrays. 

GI,  GQ  -  I  and  Q  channel  outputs  of  the  pulse  shaping 
filter. 

VlMAX,  VI  -  Variables  used  in  power  calculations. 

RI,  RQ  -  Arrays  that  store  I  and  Q  values  of  z(k) . 
RAI,  RAQ  -  Arrays  used  in  wrap  around  processing. 
PWRI,  PWRQ  -  Power  in  the  I  and  Q  components. 

FILES  USED; 


Zll.dat  -  data  file  that  stores  the  in-phase  values. 
ZQQ.dat  -  data  file  that  stores  the  quadrature  values. 


VMAX  :*  255  VLEN  :=  0  .  .VMAX  Progreun  PNCODE.MCD 


BMAX 

•  w 

•  ^ 

7 

b  ;=  0  . .BMAX 

TMAX 

•  * 

30  t  :=  0  ..TMAX 

PN  :== 

(1 

-1  -1  -1  -1  1  -1  -1 

1-11 

1-1-1111 

-1  -1-111-11 

11-1 

1  -1) 

I 

0 

0 

T  :=  0 

PULS  I 

:=  0 

0 

VLEN 

VLEN 

BIT  : 

= 

0 

0 

PULSQ 

:=  0 

0 

VLEN 

0 

<t> 

<t> 

BITS  :« 

BIT  PN 

PULSI 

•  M 

•  * 

-.038 

PULSQ 

• 

•  • 

-.058 

0 

0 

T 

:=  BITS 

b+(BMAX+l) 

t  b,  t 

PULSI 

•  “ 

-.109 

PULSQ 

.001 

1 

1 

PULSI 

•  ^ 

-.104 

PULSQ 

.133 

2 

2 

PULSI 

;  = 

.033 

PULSQ 

;  = 

.249 

3 

3 

PULSI 

:  = 

.26 

PULSQ 

.244 

4 

4 

PULSI 

;  = 

.467 

PULSQ 

;  = 

.093 

5 

5 

PUISI 

«  ^ 

.573 

PULSQ 

-.13 

6 

6 

PULSI 

•  ^ 

• 

.587 

PULSQ 

;s 

-.318 

7 

7 

PULSI 

«  w 
• 

.579 

PULSQ 

-.389 

8 

8 

98 


587 


318 


POLSI 


PUISQ  :=  - 


9  9 


PU15I  := 

.573 

PUI5Q 

10 

10 

P01£I  :  = 

.467 

PULSQ 

11 

11 

PULSI 

12 

:=  .26 

PULSQ  := 

12 

POLSI 

13 

:=  .033 

PULSQ  := 

13 

PULSI 

14 

:=  -.104 

PULSQ  := 

14 

PULSI 

15 

:=  -.109 

PULSQ  := 

15 

xfft(a)  :=  [(VMAX  +1)  ] 

icfft(a) 

Put  FIT'S 
terms  (7: 

in  Ludeman's 
286)  . 

ixfft(a)  ;= 

1 

cfft(a) 

.5 

.(VMAX  +  1) 

PI  ;=  xfft(T)  P2  ;=  xfft(PULSI) 

P3  : 

=  xf ft (PULSQ) 

•  « 

II 

(PI  P2) 

P5  ;  = 

(PI  P3) 

.13 

093 

244 

.249 

133 

001 


GI  :=  Re(ixfft(P4) ) 


GQ  :=  Re(ixfft(P5) ) 


km20c 


7 


k  0  .  .kmax 


IMAX  :=  30  ITER  ;=  0  ..IMAX 
VIMAX  :=  (IMAX  +  1)  •  (kmax  +  1)  -  1 
VI  :=  0  ..VIMAX 

RI  :=  GI 

k, ITER  k+ (kmax+1)  ITER 


RQ  :=  GQ 

k, ITER  k+ (kmax+1)  ITER 


RAI  :=  GI  RAQ  :=  GQ 

k,0  k+248  k,0  k+248 


RI  :  = 

RI  + 

RAI 

k,  0 

k,0 

k,0 

RQ  :  = 

RQ  + 

RAQ 

k,0 

k,0 

k,0 

WRITEPRN 

ZII  1 

:=  RI 

dat  J 

WRITEPRN 

'ZQQ 

;=  RQ 

dat 

Wrap-around 


s 


PWRI  := 


VIMAX 


—  S  S  h  1 

X  +  1  ^ — I  ^ — I  [.  k,ITERj 


k  ITER 


PWRQ  ;= - \  \  '  [RQ  ] 

VIMAX  +  1  ^ — I  ^ — I  [.  k,ITERj 


k  ITER 


PWRQ  = 


PWRI  = 


if 

1 
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:  REVEC.MCD 


PURPOSE;  Generates  i(k). 

SELECTED  VARIABLE  IDENTIFICATION: 

RI,  NI,  JI  -  Arrays  which  contain  the  in-phase 
components  of  1/  fi/ 

RQ,  NQ,  JQ  -  Arrays  which  contain  the  quadrature 
components  of 

RII,  Nil,  JII  -  Vectors  which  contain  the  in-phase 

components  of  l,  fi,  i. 

RQQ,  NQQ,  JQQ  -  Vectors  which  contain  the  quadrature 

components  of  2,  n, 

REVECI,  REVECQ  -  Arrays  which  contain  the  I  and  Q 

components  of  C(k) . 

REVECII,  REVECQQ  -  Vectors  which  contain  the  I  and  Q 

components  of  jr()c)  . 

ZZZI,  ZZZQ  -  Bit  synchronization  processing  vectors. 

FIX^S.  ys.sp: 

REVECI.dat  -  Stores  I  channel  components  of  r(k) . 
REVECQ.dat  -  Stores  Q  channel  components  of  r(k) . 
Note;  All  other  variables  have  been  defined  previously. 


101 


VMAX  :=  247  VLEN  :=  0  .  .VMAX  Program  REVEC.MCD 

kmcuc  :  =  7  k  :  =  0  . .  kmax 

IMAX  :=  30  ITER  :=  0  ..IMAX 


RI  ;=  READPRN 

ZII  I 

RQ  :=  READPRN 

ZQQ  ] 

dat  J 

dat  J 

NI  :=  READPRN 

'noil 

NQ  :=  READPRN 

‘noiQ  1 

dat 

dat  J 

JI  :=  READPRN 

JAMI 

JQ  :=  READPRN 

JAMQ  1 

dat 

dat  J 

<ITER> 

<ITER>  <ITER>  <ITER> 

REVECI 

=  RI 

+  NI  +  JI 

<ITER> 

<ITER>  <ITER>  <ITER> 

REVECQ 

=  RQ 

+  NQ  +  JQ 

RII 

:=  RI 

k+(kmax+l) • 

ITER 

k,ITER 

Nil 

:=  NI 

k+(kinax+l)  • 

ITER 

k,ITER 

JII 

:=  JI 

k+(kmax+l) ■ 

ITER 

k,ITER 

RQQ 

:=  RQ 

k+(kmax+l) ■ 

ITER 

k,ITER 

NQQ 

;=  NQ 

k+(kmax+l) • 

ITER 

k,ITER 

JQQ 

:=  JQ 

k+Ckmax+l)  • 

ITER 

k,ITER 
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REVECII 


:=  REVECI 


k+(]anax+l)  •  ITER  k,ITER 


REVECQQ  :=  REVECQ 

k+ (kmax+1) • ITER  k, ITER 


ZZZI 

;=  REVECII 

VLEN 

mod (VLEN+4 , VMAX+1) 

ZZZQ 

;=  REVECQQ 

VLEN 

mod ( VLEN+4 , VMAX+ 1 ) 

REVECII  :=  ZZZI 

Used  for  bit 

REVECQQ  :=  ZZZQ 

synchronization 

REVECI 

REVECII 

k, ITER 

k+ (kmax+1)  •  ITER 

REVECQ  ;= 

REVECQQ 

k,ITER 

k+ (kmax+1) • ITER 

WAME  OF  PROGRAM I  PVEC.MCD 
PURPOSE ;  Calculates  the  P-vector. 

SELECTED  VARIABLE  IDENTIFICATION: 

D  -  The  desired  signal. 

RI  -  The  in-phase  component  of  I(k) . 

RQ  -  The  quadrature  component  of  i(k) . 

R  -  The  received  array  i(k) . 

U,  V  -  Processing  arrays. 

P  -  The  P-vector. 

FILES  USED: 

PVECI.dat  -  Stores  the  in-phase  component  of  the 
P-vector. 

PVECQ.dat  -  Stores  the  quadrature  component  of  the 
P-vector. 


Note:  All  other  variables  have  been  defined  previously. 


Progreun  PVEC.MCD 

First  initialize  and  define  constants, 
kmax  :  =  7  k  :  =  0  . .  ki&ax 

IMAX  :=  30  ITER  :=  0  ..IMAX  ITERl  :=  1  ..IMAX  +  1 

D  :=  1  j 


fREVECI  1 

RQ  :=  READPRNfREVECQ  I 

[  dat  J 

dat  J 

R  :=  RI  +  RQ  j 
U  :=  R 


k  IMAX 


P  :=  D  V  31  0 


X  1  k, 


ITER 


ITER 


Used  for  troubleshooting 
purposes 


P  :  = 


0 

.17 

.5 

,849 

1 

,849 

.5 

.17 


WRITEPRNfPVECI  T  :=  Re(P) 

[  dat  J 


WRITEPRNfPVECQ  1  :=  iTn(P) 

L  dat  J 
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NAME  OF  PROGRAMS;  GRIFFPV.MCD  and  CONTPV.MCD 

PURPOSE ;  Generate  adaptation  cycles  of  the  Griffiths 

filter.  GRIFFPV.MCD  performs  the  first  "IMAX" 
iterations.  CONTPV.MCD  allows  the  user  to  perform 
as  many  adaptations  as  he  desires  in  addition  to 
those  performed  by  GRIFFPV.MCD. 

SELECTED  VARIABLE  IDENTIFICATION; 


^  -  Convergence  rate  constant. 

W  -  The  weight  vector. 

RI  -  In-phase  component  of  the  received  vector. 

RQ  ~  Quadrature  component  of  the  received  vector. 

PI  -  In-phase  component  of  the  P-vector. 

PQ  -  Quadrature  component  of  the  P-vector. 

P  -  P-vector. 

Y  -  The  actual  filter  output. 

ERR2  -  the  error  signal. 

FILES  USED; 

Wl.dat  -  Stores  in-phase  values  of  the  weight  vector 

WQ.dat  -  Stores  quadrature  values  of  the  weight 
vector 

MAGY.WKS  -  Stores  the  magnitude  output  of  the  filter 


ERR2.WKS  -  Stores  the  error  output  of  the  filter. 
Note:  All  other  variables  have  been  defined  previously. 


Program  GRIFFPV.MCD 

First  initialize  and  define  constants 


kmax  :=  7 
IMAX  ;=  30 
;=  .01 
D  :=  1  +  j 


k  : =  0  . . kmax 
ITER  :=  0  . .IMAX 


ITERl  :=  1  ..IMAX  +  1 


o 

II 

k 

-1  -1  -1  1 

-1-11-111  -1  -1 

-1  -1  -1  1 

1-1111-11-1) 

Treveci  1 

RQ  :=  READPRNfREVECQ  I 

L  dat  J 

L  dat  J 

R  :*  RI  +  RQ  j 


r  pvEci 

PQ  :=  READPRN rPVECQ  "j 

L  dat  . 

dat  J 

P  :=  PI  +  PQ  j 


The  next  equation  is  the  Griffiths  algorithm 


<ITER1> 


W 


r  <ITER1-1>'1  T  < 
[r  j  W 


<ITER1-1>  r  <ITER1-1>M  <ITER1-1> 
P  -  R  R 


<ITER1-1> 


2  +  W 


<ITER>  r  <ITER>1'''  <ITER> 
Y  :=  1_R  J  W 


<IMAX> 

Y 


ERR2 

ITER 


l[. 


<ITER> 


Re(D)  -  PN 


RelY 


<ITER>j j 


Ia(D) 


<ITQt> 

PN 


2 


<ITER> 

Y 


ERR2 

ITER 


0  ITER 

IMAX 

0 

ITER 

IMAX 

WRITEPRN 

WI  1  :  = 

dat  J 

r  <IMAX+1> 
Re[w 

WRITEPRN 

WQ  ]  :* 

dat  J 

r  <IMAX+1>1 

Iin[w  J 

1  <ITER>I 

WRITEPRN  [MAGY  "I 

:=  |y  I 

WRITEPRN rERR2  ] 

:=  ERR2 

[  ms  J 

WKS 

ITER 

! 


I 
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Program  CONTPV.MCD 

First  initialize  and  define  constants 


kmax 

:=  7 

k  :=  0  . 

.kmax 

IMAX 

;*  30 

ITER  := 

0  ..IMAX 

ITERl  :=  1  ..IMAX  +  1 

:  = 

.01 

D  :  = 

1  +  j 

W 

:=  0 

k 

PN  :  = 

=  (1  -1 

-1  -1  - 

■1  1  -1 

-11-111-1-11 

-1-1-111-1111-11  -1) 


RI  :=  READPRNFREVECI  ] 

RQ  :=  READPRNfREVECQ  I 

L  dat  J 

L  dat  J 

R  :=  RI  +  RQ  j 

PI  :»  READPRNfPVECI  "I 

PQr  :=  READPRNTpVECQ  1 

L  dat  J 

[  dat  J 

P  :=  PI  +  PQr  j 

wl  :=  READPRNfWI 

dat 


wQ  ;=  READPRNfWQ 

L  dat 


<0> 

W  :=  wl  +  wQ  j 
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The  next  equation  is  the  Griffiths  algorithm 


<ITER1-1> 

+  W 


<ITER>  r  <ITER>]''‘  <ITER> 
Y  :=  [r  J  W 


<IMAX> 

Y 


<ITER1> 


W 


<ITER1-1> 


P  -  R 


f  <ITER1-1>1 T 

Lr  J  W 


<ITER1-1> 


2  V 


ERR2 


ITER 


Re(0)  -  PN 


<ITER>  f  <ITER>1'1 

B*LV  JJ 


<ITER> 


I®(D)  -  PN 


IB 


[^<ITER>j]‘ 


<ITER>| 

f  I 

■ 

0  ITER  IMAX 

f  <IMAX-H>'1 

WRITEPRNfWI  1  :=  Re[w  J 


[fWI  1  :: 

L  dat  J 


WRITEPRNTWQ  ■)  ;=  Im 

L  dat  J 


LW  J 


<ITER> 

APPENDPRN  TMAGY  "j  :=  |Y  I  o 

WKS 


iTMAGY  I  :  = 
L  WKS  J 


APPENDPRNrERR2  I  ;=  ERR2 

L  WKS  J  ITER 


■ 

ERR2 

ITER 


0  ITER  IMAX 


<IMAX+1> 

W 
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NAME  OF  PROGRAM;  NOISE. MCD 

PURPOSE;  Generates  values  5j(k)  and  SQ(k) . 

SELECTED  VARIABLE  IDENTIFICATION; 

NPWR  -  Noise  power. 

nl  -  The  in-phase  matrix. 

nQ  -  The  quadrature  matrix. 

STD  -  The  standard  deviation. 

EILE.g_UgEO; 

noiI.dat  -  Stores  the  in-phase  values. 
noiQ.dat  -  Stores  the  quadrature  values. 

Note;  All  other  variables  have  been  defined  previously. 
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This  progr2uii  generates  the  noise  vectors  used  in 
the  main  program.  Name  of  program  "NOISE. MCD” 


NPWR  := 

1.875 

kmax 

STD  :=  ' 

Jnpwr 

IMAX 

m  :=  1 

.  .12 

k  :=  0  ..kmax  VMAX  :=  247 

ITER  :=  0  ..IMAX  VLEN  :=  0  . 


nl 

k,ITER 


(md(l)) 


6 


nQ 

k,ITER 


(md(l)) 


6 


WRITEPRNfnoil  I  :=  nl 

L  dat  , 


WRITEPRNfnoiQ  1  :=  nQ 

L  *^at 


PWRI 

PWRQ 

PWRI 

PWRQ 


1 


VMAX  +  1 


1 


VMAX  +  1 


V  V  [nl  ] 

^ — I  ^ — I  [  k,ITERj 


ITER 


k  ITER 


2 


2 


127 
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Appendix  B  :  Definitions  of  Selected  Terms 


1.  misadi ustment  -  "misadjustment  in  an  adaptive  process 
is  defined  as  the  ratio  of  the  excess  mean-square  error  to 
the  minimum  mean-square  error,  and  is  thus  a  measure  of  how 
closely  the  adaptive  process  tracks  the  true  Wiener 
solution,  that  is,  a  measure  of  the  'cost  of  adaptability’  " 
(8:110).  If  the  mean-square  value  of  the  error  signal  is 
defined  as 


e  =  E{ 

then  misadjustment  is  defined 

misadjustment  = 


e2(k)  } 

as  (10:115) 

E(  6(00)  )  - 
^min 


(B-1) 


(B-2) 


2.  cvclostationary  -  A  random  process  X(t)  is 
cyclostationary  in  the  wide  sense  if  its  mean  and 
autocorrelation  are  periodic  with  some  period  T  (13:302). 


3.  phase  distortion  -  a  linear  system  introduces  phase 
distortion  when  the  phase  response  of  the  system  is  not  a 
linear  function  of  frequency  (18:66). 


4.  roll-off  factor  -  the  roll-off  factor  is  defined  to  be 
the  excess  bandwidth  used  divided  by  the  minimum  Nyquist 
bandwidth.  The  roll-off  factor  takes  on  values  from  zero  to 
one  (9:368)  . 


5.  eraodic  random  process  -  Processes  for  which  time 
averages  and  ensemble  averages  are  interchangeable  (18:266). 
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Griffiths  algorithm  and  (2)  has  a  pseudonoise  sequence  as  an 

a 

input.  The  analysis  is  conducted  over  several  jamming 
powers,  frequencies,  and  phases.  The  Griffiths  adaptive 
matched  filter  is  shown  to  converge  for  raised  cosine  pulses 
that  experience  no  distortion,  quadratic  delay  distortion, 
and  cubic  delay  distortion.  The  Griffiths  adaptive  matched 
.filter  diverges  for  pulses  that  experience  linear  delay 
distortion  even  though  the  convergence  rate  constant  is 
within  limits.  ^Throughout  the  analysis  the  P-vector  is 
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Noise  jamming  Is  shown  to  be  more  effective  .in  the  higher 
power  ranges.  A  comparison  is  made  between  the  Griffiths 
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adaptive  matched  filter  and  an  adaptive  matched  filter  that 
uses  the  IMS  algorithm.  The  degradation  in  performance  of 
the  Griffiths  filter  compared  to  an  IMS  filter  that  uses  a 
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